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ABSTRACT 


A  new  method  of  determining  ship  model  wave  resistance  from 
the  wave  pattern  generated  during  a  test  in  a  model  tank  is  outlined. 
This  method  is  based  on  a  Fourier  analysis  of  lateral  wave-slope  data 
taken  on  a  longitudinal  cut  parallel  to  the  model  path.  Results  of 
an  exploratory  test  series  carried  out  at  Webb  Institute  are  given 
and  compared  with  previous  results  using  a  different  method.  The 
comparison  is  very  encouraging  and  indicates  that  a  new  technique 
has  been  established  which  has  many  beneficial  features,  including 
potential  application  in  full-scale  research.  Details  of  an  electrical 
wave-slope  measuring  probe  for  use  in  the  model  tank  are  included  in 
an  appendix. 
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nomenclature 

ymbols 

constants,  functions,  or  points 

wave  amplitude,  ft. 

modeL  beam,  ft. 

breadth  of  model  tank,  ft. 

Fourier  Cosine-transform;  also  2 
resistance  coefficient  =  R/^2.SVm 

block  and  prismatic  coefficients 

distance  of  cylinder  and  model  off 
tank  centerline,  ft. 

energy  in  wave  system,  ft.  lbs. 

Froude  Number  -  vm /-/I^ 

2 

acceleration  of  gravity,  ft. /sec. 
model  draft,  ft. 

depth  of  water  in  model  tank,  ft. 

,  -1 

wave  number ,  ft. 
constant 

model  length,  ft. 

2 

pressure,  lbs. /ft. 
resistance,  lb. 

element  of  length  along  path  AB,  ft. 
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Fourier  Sine^transform;  also  wetted 
surface,  ft. 

fluid  velocity  components  in  the 
x,  y,  z  directions,  ft./£ec. 

Eggers'  parameters,  see  Equation  (1$ 

model  speed,  ft. /sec. 

fluid  velocity  normal  to  AB,  ft. /sec. 

o 

model  weight,  lb. 

energy  transfer  into  wave  system  from 
rear  plane,  ft.  lbs. /sec. 


Cartesian  co-ordinates,  ft.  (see 
Figure  2)  is  vertically  upwards. 

angle  of  path  with  the  x-axis,  radians 

2  2 

speed  constant  gL/Vm  =  1/Ff 
displacement,  lbs. 
wave  elevation,  ft. 


radians 


lateral  wave  slope 
wave  direction  angle 
wave  length,  ft. 
index  -  1,  2,  •  •  •  •  7/mtUi 

3 

density  of  fluid,  slugs/ft. 

velocity  potential ,  sujh  that 
u  ■  ”  °  »  etc,»  ft*  /sec* 

wave  coordinate  ■  x  cos0  +  y  sin0,  ft. 


ohms  electrical  resistance 
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I.T.T.C. 
t,  f,  r,  w 


centerline  wave  properties 
lateral  slope 

International  Towing  Tank  Conference 

total,  frictional,  residual,  and 
wave  parts 

(Illustrated  for  "Q") 

time  rate  of  change  of  Q 

Q  per  unit  distance  in  the  x  direction. 
Primes  are  also  used  to  denote  variations 
in  a  quantity. 

average  value  of  Q  (see  footnote  page  6  ) 
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INTRODUCTION 


In  recent  years  several  methods  have  been  developed  to 

determine  ship  model  wave  resistance  experimentally  in  the  model  tank 

from  measurements  of  the  wave  system  generated  by  the  model  during  a 

test.  Interest  in  this  subject  was  stimulated  by  a  conference  at  Ann 
1 

Arbor  sponsored  by  the  Office  of  Naval  Research.  The  methods,  several 
of  which  have  been  found  to  be  successful,  stem  from  essentially  the 
same  assumptions  of  linearized  potential  wave  theory,  but  differ  in 
the  type-'of  wave  characteristic  measured  (e.g.  wave  elevation,  slope, 
or  force  on  a  fixed  object)  and  the  type  of  path  along  which  data  are 
taken  (e.g.  a  longitudinal  or  transverse  cut,  relative  to  the  model 
path).  Most  depend  upon  a  Fourier  analysis  of  the  data  and  yield  the 
spectrum  of  the  wave  pattern  as  well  as  the  wave  resistance  as  the 
result.  It  should  be  noted  that  it  is  not  necessary  to  assume  the 
validity  of  any  ship-wave  generation  theory  (e.g.  thin  ship  theory) 
for  the  derivation  of  these  methods.  An  up-to-date  picture  of  this 
work  including  a  complete  bibliography  of  papers  related  to  this 
subject  is  contained  in  Reference  (2). 

The  work  done  by  the  author  at  Webb  Institute  over  the  past 
3  4 

few  years  *  has  primarily  concerned  his  so-called  "X-Y  Method"  which 
uses  the  x  and  y  components  of  the  lateral  force  exerted  by  the  wave 
system  on  a  long  fixed  vertical  circular  cylinder,  and  which  does  not 
require  a  Fourier  analysis  but  yields  the  wave  resistance  directly. 

A  major  advantage  of  this  scheme,  as  compared  to  those  using  wave 
elevation  data,  was  found  to  be  the  finite  length  or  "self  truncation" 
of  the  record.  This  results  from  the  rapid  decay  of  the  y  force  signal 
after  passage  of  the  model,  before  wave-reflection  from  the  side  walls 
of  the  tank  becomes  important.  Recently  a  method  with  a  similar 

Raised  numbers  denote  References  listed  beginning  on  page  40. 
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advantage  has  been  tried  utilizing  lateral  wave-slope  data  taken  along 

a  longitudinal  cut.  This  method  was  first  suggested  along  with  a  number 

5 

of  others  by  Sharma  in  1963  and  yields  the  wave-slope  spectrum  as  well 

as  the  wave  resistance.  Tests  run  at  Webb  Institute  in  June  1965,  just 

before  the  author  left  for  a  one  year  National  Science  Foundation 

Fellowship  in  Hamburg,  Germany,  have  now  been  analyzed,  and  this  method 

has  been  found  to  be  quite  convenient  and  successful  in  the  smaller 

facility  at  Webb.  An  attempt  to  include  this  method  among  others  during 

2 

a  test  series  in  the  Hamburg  model  tank  did  not  meet  with  similar  success, 

mainly  because  of  the  difficulties  encountered  there  in  connection  with 

the  slope  probe  instrumentation  problem.  However,  in  connection  with 

2 

the  latter  effort  computational  and  theoretical  evidence  was  found 
which  indicated  this  method  to  be  consistent  with  some  of  the  others 
which  were  used  in  the  actual  tests.  This  paper  will  outline  the  lateral 
wave-slope  method  and  underlying  theory  briefly  and  present  the  test 
results  and  experience  gained  to  date  in  its  use. 


THEORY 


The  basic  underlying  theory  of  wave  survey  methods  in  general 

can  be  found  in  several  references,  the  latest  being  that  by  Eggers, 

2 

Sharma,  and  Ward  .  In  this  section  a  brief  outline  will  be  given  of  the 
theory  needed  for  the  lateral  wave-slope  method  in  particular,  with  emphasis 
on  understanding  the  significance  of  the  various  steps  involved  in  carrying 
it  out. 

Momentum  Analysis 

ic 

We  consider  the  situation  shown  in  Figure  1  of  a  model  traveling 

at  constant  speed  V  along  the  centerline  of  a  long  tank  of  breadth  b  and 

m 

depth  h,  and  a  set  of  axes  x,  y,  z  moving  with  the  model.  We  consider 

potential  flow  and  small  surface  disturbances,  •  Consider  a 

volumetric  region  I,  fixed  in  space,  whose  side  boundaries  are  formed  by 

the  lines  ABCDC'B'  projected  down  from  the  free  surface  except  in  the 

region  of  the  model  sides  (where  it  is  formed  by  them)  and  whose  bottom 

surface  is  the  tank  bottom,  both  of  the  latter  within  the  lines  ABCDC'B'. 

We  consider  path  with  symmetry  about  the  tank  centerline,  and  therefore 

only  the  region  ABCD  on  one  side  of  the  tank  need  be  considered  and  the 

results  doubled.  The  line  CD  is  drawn  sufficiently  far  ahead  of  the  model 

so  that  there  is  no  disturbance  there.  The  line  AB  in  general  forms  a  curve 

as  shown.  The  model  moves  ahead  at  the  speed  V  and  creates  waves  within 

m 

region  I.  This  wave  system  is  often  idealized  to  the  concept  of  a  transverse 
system  0  and  a  divergent  system  @  contained  within  the  cusp  lines  as 


Figures  are  numbered  in  order  and  located  at  appropriate  places  in 
the  text.  An  index  of  Figures  and  Tables  is  given  on  page  iii. 
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shown  in  Figure  1  and  is  often  referred  to  as  the  Kelvin  wave  pattern. 
The  wave  system  grows  within  this  region  and  therefore  the  energy  in 
the  wave  system  increases  at  the  rate: 


E 


E  VM 


(i) 


where  E'  is  twice  the  average  total  wave  energy  s?  a  strip  of  unit 
width  in  the  x  direction  along  the  line  AB  extended  av,:*n  through  the 
depth  h.  This  growth  of  energy  must  be  supplied  by  the  work  done  at 
the  boundaries  of  the  region,  this  rate  of  doing  work  being  given  by 
integrals  over  the  boundaries  of  the  fluid  of  the  pressure  p  times 
the  velocity,  Vn,  normal  to  the  boundary.  There  are  no  contributions 
to  W  from  the  tank  sides  BC  and  B'C'  or  the  front  plane  CDC'  or  the 
tank  bottom  due  to  the  normal  velocity  being  zero  in  each  case,  nor 
from  the  free  surface,  due  to  the  pressure  there  being  zero  ambient. 
Along  the  model  sides  there  is  work  being  put  in,  due  to  the  action 
of  the  pressures  there  in  conjunction  with  the  forward  motion  of  the 
model,  at  the  rate  R  V  ,  where  R  is  the  wave  resistance,  which  is 
the  only  resistance  possible  in  the  present  case  in  which  an  ideal 
fluid  and  small  disturbances  are  assumed.  Along  the  prismatic 
surface  AB  we  have  the  work  done  by  the  waves  in  region  II  on  region 
I,  i.e.  the  "work  done  on  the  rear  plane": 

r  r? 


P  Vh  dl 


-h 


(2) 


where  p  is  the  pressure,  is  the  fluid  velocity  normal  to  the  surface 
AB,  and  ds  is  the  element  of  distance  along  AB.  Due  to  symmetry  we 
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•  • 

have  WQ'A  «  W*e  •  For  the  general  case  shown  where  the  curve  AB 
makes  an  angle  (  -  o<  )  locally  with  the  x  axis  we  have: 


where :  i*  =  fan  o<  and, 

da  ' 

u  and  w  are  the  fluid  velocities  in  the  x  and  y  directions. 
The  pressure  p  is  given  by  the  linearized  Bernoulli  relationship: 


p  -  p.  +  r uV*  -  i°9 1 


(4) 


where  terms  of  the  order  of  the  fluid  velocities  squared  have  been 
neglected  in  relationship  to  the  unsteady  pressure  term  P  U  . 
Since  the  value  of  W  for  use  in  the  energy  relationship  must  be  an 
average  value  with  respect  to  time  at  a  fixed  plane,  the  pressure 
times  velocity  terms  represented  by  the  first  and  third  terms  in 
Equation  (6)  will  average  to  zero.  Thus  we  have: 


r 


u(u -vfttnc*)  d  Z 


Since  this  is  the  only  work  term  other  than  that  due  to  the 
force  we  have,  equating  the  total  work  put  into  region  1  to 
average  rate  of  growth  of  energy  in  the  wave  system  in  this 


(5) 

towing 

the 

region: 


wAa  =  e'vm 

Rw  =  E  -  %/Vm 


(6) 


*  The  time  average  Q  of  a  quantity  Q  is  defined  as: 


3  «2T 


Qtt)  dt 
4 


where  At  is  the  period  of  fluctuation  of  Q  or  the  limit  can  be  taken 
as  A"t  -»  ®o 

-6- 


The  total  energy  E'  Is  equal  to  the  sum  of  the  potential  and  kinetic 
energies  per  unit  distance  in  the  x  direction: 


E  = 


[  f  d!j  +  fl dy  [  (u2  +  v*+wz)  dz 

Mb  Mb  J-t, 


(7) 


where  fs  is  the  wave  elevation  and  w  is  the  fluid  velocity  in 

the  2  direction.  Combining  Equations  (5), (6),  and  (7)  gives  the 
"1st  Havelock  Formula": 


Rw  =  fg  f f2 d<j  +  e\ f(ua+v2+w2)  d 

JAB  J  MB  J-h 

- 2. /»  Ju(u-v£an<*)  di 


(8a) 


for  the  general  path  AB.  Two  special  cases  are  of  interest:  For 
plane  perpendicular  to  the  direction  of  motion  we  set  ©<  «  0  t 
giving  i 

•V*  rt 

(.^+yz+WZ)  dz  (8b) 

-o  '  'J  'J 


f  a  1*5  « 

r‘/i 

Rw  •  rsj 


o 


This  is  the  form  used  for  a  transverse  cut  method.  For  a  longitudinal 
cut  method  we  set  C*  ,  giving: 


R 


w 


2y°  c/k  Ju  V  d  2 


(8c) 
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Fourier  Wave-Spectrum  Model 


A  study  of  Equation  (8a)  will  convince  the  reader  that  this 
relationship  does  not  form  a  convenient  basis  for  an  experimental 
method,  as  it  requires  determination  not  only  of  the  wave  elevation 
along  the  path  AB  but  of  the  fluid  velocities  u,  v,  w  throughout  the 
depth  of  khe  fluid  along  the  vertical  plane  formed  by  the  line  AB 
extended  in  depth  into  the  fluid.  We  wish  to  use  potential  wave 
theory,  assuming  free  gravity  waves  of  small  amplitude,  in  effect  to 
predict  these  fluid  velocities  from  knowledge  of  the  wave  elevation 
(or  in  the  present  case  from  knowledge  of  the  wave  slope  in  the  y 
direction  T  y  This  will  allow  us  to  "trade-off"  more  detailed 

knowledge  of  the  function  ^  on  the  surface  for  the  next-to-impossible 
integration  in  depth.  This  would  be  quite  simple  for  a  single  component 
wave  of  amplitude  A  and  length  X  traveling  in  a  direction  0  to 


the  x  axis: 


I  =  A  cos  kco 


where :  <0-  XCOS©+^Sin6 

k  *  * 

Potential  gravity  wave  theory  would  then  predict: 


u  = 

V  = 

W  -  - 


where: 


••III 


6  cos  kco 


is  the  wave  speed  in  deep  water  which  is  assumed.  The 


tt 
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wave  system  which  is  likely  to  be  generated  by  even  a  simple  traveling 
disturbance  is  quite  complex,  however^.  For  this  reason  we  must  use  a 
spectrum  representation: 


r% 


t 


+■  B<@)  sin  kw  de  (u) 


zy2 


This  can  be  seen  to  be  the  superposition  of  a  large  number  of  waves  of 
amplitude  A(9)d0  traveling  in  different  directions  with  a  varying  phase 
relationship  allowed  by  B(9).  A  similar  spectrum  form  exists  for  <|> 
and  therefore  the  velocities  u,  v,  w  are  known  in  terms  of  the  spectrum 
functions  A  and  B.  Moreover  simple  differentiation  would  result  in 
the  or  "lateral  wave-slope"  spectrum: 

r  % 

[AijCS)  Cos  kco  +-  BjjtSl  Si  n  kco]  dS  <12) 

where  A*9)  =  ksi*©  BO)  ,  etc. 

It  can  therefore  be  seen  that  knowledge  of  the  spectrum  of  one  characteristic 
of  the  wave  system  (say  )  will  in  theory  allow  any  other  to  be  derived 
analytically.  This  is  an  important  point  which  will  allow  one  to  pick 
the  most  convenient  data  to  measure  experimentally,  regardless  of  which 
type  of  spectrum  is  needed.  Some  typical  wave  spectra  can  be  seen  in 
Figures  5(a)-(e),  and  further  discussion  on  the  matter  of  wave  patterns 
can  be  found  in  Reference  6. 
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Stationary  Wave-Pattern  Requirement 

Another  very  Important  fact  about  the  ship  wave  system  Is  that 
It  must  be  stationary  with  relation  to  the  ship.  The  corresponding 
relationship  between  R  and  0,  which  can  be  derived  either  geometrically 
or  from  mathematical  reasoning  ,  is: 


k  m  k.  secz6 

where:  k„  -  8/ 

Other  useful  relationships  are: 

>  «  X0  Cos20 

► 

C  -a  C0  Co  S  B 


(13) 


(1A) 


where  c  is  the  wave  speed,  and  the  zero  subscript  denotes  the  wave  component 

for  9  ■  0,  or  the  "centerline"  wave.  It  should  be  noted  that  all  other 

waves  are  shorter  and  slower  than  the  centerline  wave,  and  that  c-  ■  V  . 

’  0  m 


Eggers  wy,  uv  Parameters 


Using  the  foregoing  relationships  the  argument  (k«)  of 
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the  wave  functions  can  be  written: 

kco  —  Wy  X  +*  U  v  y 

where:  Wv  =  k0  Sec© 

«  koSec0fan@ 

The  factors  w-v  and  u^are  the  convenient  "wave  numbers  in  the 

y  8 

x  and  y  directions  introduced  by  Eggers  ,  and  are  used  as  parameters 
in  the  analysis  to  follow.  For  this  purpose  however  they  are  non- 
dimensionalized  as  follows: 

W y  -  Wv/ 

Uy  -  Uv/k. 

This  form  will  be  used  with  the  primes  not  written  in  this  text.  It 
is  seen  that  row  there  is  only  one  parameter  in  the  spectrum  representation; 
either  9,  wv  or  u^with  the  other  two  derived  from  it.  Many  different 
forms  of  the  previous  equations  and  those  to  follow  are  therefore  possible. 

Wave  Resistance  Formula 

With  the  wave  spectrum  model  now  established,  the  1st  Havelock 
Forumla  (8)  can  now  be  converted  to  integrations  involving  only  the  function 
that  is  the  surface  wave  elevation.  The  resulting  equations  will  appear 
either  (a)  in  terms  of  the  wave  spectrum  functions  A  and  B  or  (b)  in  terms 
of  the  related  Fourier  transform  coefficients  (denoted  S,  C,  etc.)  of  aaae 
type  of  wave  data  (  £  , J  etc.)  along  some  path  (longitudinal  or  transverse 
cut,  etc.),  or  (c),in  the  special  case  of  the  X-Y  Method,  in  terms  of  a 


direct  integral  of  a  ^-function  along  a  path.  The  latter  is  fully 
discussed  in  Reference  (3).  The  first  mentioned  is  the  classical 
•'2nd  Havelock  Formula"9: 


(a2+b2)  cos3e 


de 


(17) 


which  is  extremely  simple  in  form  and  requires  only  the  magnitude  of 
the  wave  spectrum  A(9),  B(0)  to  be  known.  The  second-mentioned  is  very 
closely  related  to  it  and  gives,  for  lateral-slope  data  taken  on  an 
infinite  longitudinal  cut  : 


where : 


dUv 

V(2.w/-l) 


1 


S* 

c, 


* 


dx 


(18) 


(19) 


are  the  sine  and  cosine  transforms  of  the  lateral  slope  data  along  a 
line  y  *  constant.  The  infinite  limits  are  in  practice  replaced  by 
a  finite  integration  over  the  extent  of  the  signal  which  for  the 
lateral  slope  data  is  finite  due  to  its  self • truncation  characteristics 
noted  previously.  It  is  interesting  to  note  that  the  6  form  of  (19): 


R 


N N 


r% 

(s*  4-  Cj2)  co s3e  de 

o 


(20) 
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is  very  similar  in  form  to  (17).  The  Eggers'  parameters  wv 

are  more  convenient  to  use  however  in  the  computations. 

Choice  of  Path 


It  should  be  noted  that  the  results  R  and  S  ,  C  should  r.ot 

w  y  y 

depend  on  the  path  chosen;  however  both  theoretical  reasoning  and 

computational  experience  show  that  the  line  ym  constant  should  be 

neither  too  near  the  model  or  the  tank  wall  for  accurate  results.  The 

latter  is  due  to  tank  wall  reflections,  not  considered  in  the  foregoing, 

and  the  former  due  to  the  presence  near  the  model  of  a  "local"  wave 

system  for  which  the  wave  spectrum  model  as  given  in  Equations  (11)  and 

(12)  is  not  a  valid  representation.  Further  guidance  on  this  question 

can  be  obtained  in  Reference  (2).  In  the  case  of  lateral  slope  data, 

for  reasonable  distances  of  the  probe  from  the  model,  tank  wall 

reflection  in  a  tank  of  normal  breadth  does  not  occur  until  the  record 

has  truncated  itself;  an  advantage  which  this  method  shares  with  the 
3 

X-Y  Method  ,  and  which  is  not  in  general  true  of  other  longitudinal 
2 

cut  methods  .  If  the  tank  breadth  is  a  critical  factor,  some  advantage 
can  be  gained  by  employing  "off-centerline  towing"  as  described  in 
Reference  (3)  and  shown  in  Figure  2.  A  longitudinal  cut  is  accomplished 
most  easily  by  employing  a  fixed  measuring  location  over  the  water  surface 
near  the  tank  centerline  and  letting  the  model  pass  by. 

In  the  foregoing,  port-starboard  symmetry  of  the  wave  system 
and  therefore  of  the  spectrum  has  been  assumed  and  therefore  only  one 
path  on  one  side  of  the  model  is  necessary.  If  such  symmetry  cannot 
be  assumed  (e.g.  a  heeled-yacht  test),  two  paths  will  be  necessary  in 
practice  to  cross  the  wave-pattern  in  the  "Kelvin  sensed  thus  re¬ 
quiring  simultaneous  recording  from  two  probes. 
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EXPERIMENTS 


Tests 


Experiments  were  carried  out  on  June  4,  1965  at  Webb  Institute 
on  a  5 ' -0"  model  of  the  Series  60,  0.60  Block  Coefficient  form.  Particulars 
of  this  hull  form  are  given  in  Table  I. 


Table  I:  Characteristics  of  the  Series  60  Model 
(See  Reference  10) 


Length  Between  Perpendiculars  L 

Waterline  Length 

Beam  B 

Draft  H 


Displacement 
Wetted  Surface 
Block  Coefficient 
Prismatic  Coefficient 


A 


s 

c. 


Displacement-Length  Ratio 


¥z240{V,cc)S 


Stimulators:  1/8"  Diameter  x  .050"  height  pins  spaced 
0.275"  placed  4"  aft  of  bow 


5.000  ft. 
5.084  ft. 
0.667  ft. 
0.267  ft. 
31.170  lbs. 
4.262  ft.2 
0.600 
0.614 

122 


FW 


The  model  basin  at  Webb  Institute  is  93  ft.  long  with  a  10  ft.  wide  by 
5  ft.  deep  rectangular  section  and  a  carriage  on  an  overhead  rail.  Tht 
model  was  towed  at  a  distance  =  6in.  off  the  tank  centerline  with 
a  slope  measuring  probe  mounted  at  a  distance  6-7/8  in.  off  the 
tank  centerline  on  a  fixed  arm  at  about  the  mid-length  of  the  tank  but 
on  the  opposite  side,  as  shown  in  Figure  2.  This  maximizes  the  length 


of  the  resulting  record  along  the  line  AB  before  wave  reflection  from 
the  sides  of  the  tank  take  place. 

Records  were  taken  of  the  lateral  slope  %  of  the  waves  caused 
by  the  model  measured  by  a  three-wire  resistance  type  slope  gage  developed 
in  conjunction  with  the  work  and  described  in  detail  in  Appendix  A.  This 
instrument  was  designed  to  produce  a  linear  reading,  on  a  two  channel 
Brush  RD-4622-00  amplifier  and  Mark  II  recorder,  proportional  to  the  wave- 
slope  in  the  y  direction  over  a  range  of  -  0.20  radians  and  to  ignore 
changes  in  the  wave  elevation  over  the  maximum  range  encountered.  Three 
repeat  runs  each  were  made  at  three  different  sf  »eds  over  the  range 
covered  pieviously  for  this  model  in  the  "X-Y*'Method"  tests,  in  order 
to  provide  an  exploratory  series  to  examine  the  method  itself.  These 
conditions  as  well  as  the  resulting  resistance  values  are  shown  in  Table  II. 
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Run  on  6/4/65  on  the  Series  60,  0.60  Block  Model 


Run  No. 

Model  Speed 

Froude 

Number 

Effective 
Tank  Width 

No.  Calc. 
Steps 

Resistance 
Coeff t. 

Vm,  ft. /sec. 

F 

r 

b  „  ft. 
ef  f  * 

V  max 

3 

10  C 

w 

D-l 

3.54 

.279 

12.2 

51 

.655 

D-l 

3.54 

.279 

12.2 

76 

.666  =  +27. 

D-2 

3.51 

.277 

12.0 

51 

.  9 

D-3 

3.50 

.276 

12.0 

51 

.494 

G-l 

4.29 

.338 

18.0 

51 

1.576 

G-2 

4.30 

.339 

18.1 

51 

1.480 

G-3 

4.30 

.339 

18.1 

51 

1.533 

J-l 

5.09 

.401 

25.3 

51 

5.759 

J-2 

5.07 

.400 

25.1 

51 

5.772 

J-3 

5.09 

.401 

25.3 

51 

5.994 
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A  sample  record,  Run  G-2,  is  shown  in  Figure  3.  It  can  be  seen 
that  this  record  is  free  from  electrical  or  mechanical  "noise".  The 
event  mark  indicates  passage  of  the  model  bow  past  the  measuring  station. 

The  recorder  was  run  at  its  maximum  paper  speed  of  125mm/sec.  to  obtain 
an  easily  readable  record.  The  calibration  factor  for  the  tests,  converted 
to  the  sensitivity  'jetting  of  this  tape,  was  85.0  run  radian,  corresponding 
to  a  full  scale  tape  reading  of  -  0.17  radians.  Less  than  1.0  mm  deflection 
was  noted  over  a  change  of  elevation  of  *  1.0".  As  can  be  seen  in  Figure  3, 
the  record  was  practically  self-ti uncating  before  wave  reflection  from  the 
tank  walls  took  place.  The  tests  v.-'re  all  accomplished,  including  set-up 
of  equipment  and  calibrations  before  and  after,  in  one  working  day. 

Although  turbulence  stimulation  is  not  as  important  in  this  type 
of  testing  as  in  the  case  of  ordinary  resistance  tests,  a  turbulence 
stimulation  pin  configuration  found  suitable  for  this  model  in  the  high 
speed  range  was  fitted  as  shown  in  Table  I,  however  the  tank  water  was 
not  heated  but  was  at  the  ambient  temperature  of  approximately  72°F.  The 
regular  side  cylindrical  wave  dampers  and  end  corner  beaches  were  in  place 
during  the  tests.  Speed  measurement  was  accomplished  by  means  of  a  timer 
acting  between  switches  on  the  rail  35.0  ft.  apart,  and  the  model  was 
towed  by  an  arm  on  a  fixed  strut  on  the  carriage. 

Analysis 

Prior  to  actual  analysis,  each  tape  (e.g.Fig.  3)  was  divided  into 
2.5  mm  spacing  the  divisions  shown)  and  numbered  and  the  data  read  off 
and  transferred  to  IBM  cards,  thus  providing  the  raw  data  in  a  usable 
"digitized"  form.  While  this  was  of  necessity  done  by  hand,  the  instrumentation 
would  be  easily  adapted  to  automatic  digitizing  by  suitably  designed  tape 


I 


P  »5UR€  3’ 

Rccdrd  (Q-2) 


apparatus  with  a  "playback"  feature  to  realize  the  scanning  rate  of  50/sec. 
corresponding  to  the  hand  calculation. 

The  analysis  itself  consists  of  two  steps.  First  is  the  carrying 
out  of  the  two  x  integrations  over  each  record,  given  by  Equation  (19), 
at  a  series  of  values  of  the  "x-wave  number"  wy  to  give  the  Fourier  sine 
and  cosine  transforms  and  as  a  function  of  wv  .  These  are  in 
themselves  a  result  and  can  be  plotted.  The  second  step  is  the  integration 
of  the  sum  of  the  squares  of  these  with  the  weighting  function: 


K  Cwv)  = 


1 

w/(zw/-l  ) 


(21) 


in  accordance  with  Equation  (18),  to  give  the  wave  resistance  R  for 

w 

the  tape  analyzed  at  the  speed  of  the  run.  These  results  are  non- 
dimensionalized  in  terms  of  the  resistance  coefficient  C  and  the  Froude 
number  Ff  given  by: 


C  =  R/>/i  s  vjf 

K  - 


j 


(22) 


where  S  is  the  wetted  surface  and  L  the  length  of  the  model.  The 
integrations  are  done  by  trapezoidal  rule  with  steps  of  u^given  by: 


'Jv  * 


1?  V 


(23) 


where  is  a  hypothetical  effective  tank  width  corresponding  to  the 
exact  Eggers  series-spectrum  model  of  the  same  spacing.  In  the  present 
case  a  choice  of: 


K  b<ff  =  i  o  it 
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(24) 


was  made,  thus  giving  steps  A  uv  of  0.10.-,  The  corresponding  steps 
of  are  not  even  and  are  given  by: 


W y  -  \  (  I  +  )  <25> 

The  Fortran  computer  program  designed  to  carry  out  the  above 
calculation  is  described  in  Appendix  B.  A  cut-off  corresponding  to  51 

was  found  to  be  sufficient  in  terms  of  no  further  contribution  to  the 
resistance  values  beyond  this  point,  although  the  spectrum  itself  continues, 
because  of  the  strong  effect  of  the  weighting  function  K(w^)  in  Equation  (18); 
as  shown  in  Table  II  extension  to  "]/  =76  in  the  case  of  tape  D-l  gave  only 

a  27.  increase  in  the  result.  The  calculations  were  run  on  an  IBM-1130 
computer.  Analysis  time  for  the  final  set  of  calculations  was  about  3 
minutes  per  ship  run. 

Results 


The  resulting  wave  resistance  coefficients  from  this  set  of  tests 

3 

are  given  in  Table  II  and  compared  with  the  previous  results  using  the 
X-Y  Method,  which  are  shown  as  a  faired  line,  in  Figure  4.  Also  included 
in  this  figure  is  the  residual  resistance  coefficient  Cr  derived  from  a 
previous  total  resistance  test  of  the  same  model  at  Webb.  It  can  be 
seen  that  excellent  agreement  is  present  between  the  lateral  slope  results 
and  the  faired  line  at  the  lowest  and  middle  speed  tested,  and  that  the 
lateral  slope  results  are  higher  than  the  XY  at  the  highest  speed,  but 
still  less  than  the  residual.  This  latter  difference  is  considered  to  be 
in  the  right  direction,  as  the  XY  results  are  known  to  be  in  error  above 
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a  Froude  number  of  0.34  for  these  tests  due  to  the  XY  cylinder  being 

of  too  small  a  diameter.  Moreover,  it  can  be  argued  that  the  measured 

wave  resistance  should  approach  the  residual  in  the  range  of  high 

Froude  numbers.  The  basic  validity  of  the  X-Y  Method  has  been  further 

2 

established  by  means  of  recent  tests  in  Hamburg,  Germany  on  a  larger 
model.  Thus  from  the  limited  number  of  results  obtained  so  far  the 
lateral  slope  method  appears  to  be  a  valid  technique  for  measurement 
of  wave  resistance  in  the  model  tank,  giving  results  as  good  as  and 
perhaps  better  than  the  other  methods  used  to  date. 

Several  of  the  lateral  wave  slope  spectrum  results  are  shown 
in  Figures  5(a)-(e).  While  these  cannot  be  compared  with  any  other 
established  result,  certain  known  features  of  such  spectra  can  be 
examined.  One  is  the  presence  of  zeros  in  the  ideal  fluid  case  or 
near-zeros  in  the  actual  case  of  a  real  fluid  of  small  viscosity. 
Physical  reasoning  shows**  that  these  should  be  spaced,  when  plotted 
against  a  base  of  wv  ,  approximately  at  an  interval  of: 


A  Wv  = 


(26) 


where Xo=  It  can  be  seen  1°  Figure  5(b)  for  the  case  of  G-2 

that  this  is  approximately  true.  The  trend  of  movement  of  these  "zeros" 
away  from  the  axis  as  the  speeds  are  increased  might  be  the  result  of 
worsening  separation  at  the  higher  speeds.  In  addition  the  consistency 
of  the  resulting  spectra  can  be  noted  for  the  highest  speed  tested  in 
Figures  5(c)-(e);  this  is  not  too  bad  considering  the  accuracy  of 
measurement  possible  and  the  complexity  of  the  analysis.  An  additional 
factor  which  might  contribute  to  such  variations  is  a  port-starboard 
fluctuation  of  the  wave  system  which  has  sometimes  been  noticed  during 
model  tests. 

The  weighting  function  K(wy  )  in  Equation  (11)  is  also  graphed 
in  Figure  5(b).  From  this  we  can  see  that  only  a  portion  of  the  spectrum 
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that  for  0<6<601>  approximately,  seems  to  be  needed  for  the  accurate 
determination  of  the  wave  resistance.  This  "paradox"  is  related  to  the 
fact  that  each  part  of  the  ship  wave  spectrum  is  analytically  related 
to  each  other  as  discussed  in  Reference  (6). 

It  should  be  mentioned  that  the  above  results  were  preceeded 

by  an  independent  calculation  by  S.  D.  Sharma  at  the  Hamburgische 

Schif fbau-Versuchsanstalt  in  Hamburg,  Germany  on  the  runs  D-l  through 

D-3  which  yielded  results  pratically  identical  with  those  exhibited  here. 

An  interesting  feature  in  one  of  these  results,  which  was  carried  out 

to  a  much  larger  range  of  wv  than  necessary  to  calculate  the  wave  resistance, 

is  the  persisting  nature  of  the  oscillations  in  the  spectrum  at  a  constant 

interval  close  to  the  value  Y  • 

'  °o 


-29- 


CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  comparison  of  the  results  of  the  exploratory  tests 
described  herein  with  other  related  results  in  the  area  of  wave  survey 
testing,  it  is  concluded  that  the  lateral  wave -slope  method  is  a  valid 
experimental  technique  for  determining  the  wave  resistance  and  corresponding 
wave  spectrum  of  a  ship  model  in  a  typical  model  tank.  Both  the  necessary 
instrumentation  and  the  tests  and  analysis  were  accomplished  with  reason¬ 
able  ease  and  accuracy  in  the  relatively  small  facility  at  Webb  Institute. 
Further  investigations  of  the  method  and  technique  should  include  dynamic 
calibration  of  the  slope  probe  as  well  as  experiments  at  varying  lateral 
positions  of^ the  measuring  station.  Experience  with  other  longitudinal 
cut  methods  has  shown  that  the  latter  might  have  an  influence  on  the  results. 
In  addition  it  would  be  useful  to  apply  the  method  to  other  types  of  ship 
forms  and  other  vehicle  types  to  gain  further  experience  with  these  and 
with  the  measuring  problems  involved. 

It  is  considered  that  the  method  can  now  be  applied  to  solving 
useful  and  practical  problems  for  which  it  is  ideally  suited  as  well,  such 
as  optimizing  bulbs  and  other  hull  form  variations.  Provision  of  this 
equipment  can  "stretch"  the  capability  of  a  model  tank,  especially  one 
in  the  smaller  size  range.  Also,  since  this  method  has  potential 
application  in  neai  or  full-scale  testing,  some  development  should  be 
considered  in  this  direction.  Such  studies  should  include  consideration 
of  ambient  sea  noise  as  well  as  more  practical  matters  such  as  wave- 
slope  measuring  buoys,  etc. 

Although  the  measurements  and  analysis  in  the  present  series  are 
not  considered  burdensome,  for  more  extensive  use  of  this  technique  on 
a  regular  basis  automatic  digitizing  including  print  or  punch-out  of  the 
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recorded  slope  signal  would  be  recommended.  This  would  allow  more  runs 
to  be  made  and  more  different  conditions  analyzed  as  well  as  producing 
almost  immediate  results  after  a  test. 
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APPENDIX  A:  DEVELOPMENT  OF  THE  WAVE-SLOPE 
MEASURING  PROBE 


A  successful  wave-slope  measuring  probe  based  on  the  resistance 

wire  principle  was  developed  in  1965  for  use  at  Webb  in  conjunction  with 

the  2KC  excitation  from  a  two-channt 1  Brush  RD-4622-00  amplifier  and 

Mark  II  recorder.  The  probe  and  the  calibration  scheme  is  shown  in 

Figure  A-l  and  the  corresponding  circuit  diagram  is  shown  in  Figure  A-2. 

The  scheme  is  an  adaptation  of  the  wave  height  measuring  probe  previously 

3 

developed  in  connection  with  this  effort  ,  and  employs  three  parallel 
stainless-steel  wires  of  .004  inch  diameter,  spaced  at  3/4  inch  intervals 
between  carefully  drilled  holes  in  Plexiglass  blocks  held  by  a  C-shaped 
frame.  The  frame  is  mounted  on  a  slider  which  allows  known  variations 
of  vertical  height  of  -  3"  at  approximately  V'  intervals  and  this  in 
turn  is  on  a  plate  which  can  be  rotated  in  either  direction  to  produce 
a  slope  variation  of  *  .20  radians  in  steps  of  0.10  for  the  slop'' 
calibration.  The  height  variations  are  used  to  calibrate  wave  height 
probes  mounted  on  the  same  frame  or  in  the  present  case  to  check  the 
insensitivity  of  the  wave-slope  probe  to  changes  in  wave  elevation. 

As  shown  in  Figure  A-2,  the  measuring  circuit  consists  of  a 
Wheatstone  bridge  with  two  of  the  resistances  R^  and  R^  replaced  by  the 
equivalent  path  resistance  of  the  water  between  each  set  of  wires  AB  and 
AC.  This  scheme  is  known  in  the  case  of  the  height  probe  to  give  a  good 
linear  response  to  changes  in  wave  height  over  a  range  of  several  inches 
when  the  ratio  cf  resistance: 

R3/R4  =  R2/Rl  =  K  (A_1) 

is  properly  chosen;  in  this  case  the  value  of  K  is  approximately  20  based 
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on  the  equivalent  D.C.  resistance  of  the  water  elements.  In  the  case 
of  the  slope  probe  the  elements  are  used  differentially,  as  shown.  This 
circuit  is  more  symmetric  and  actually  easier  to  balance  than  the  height- 
probe  circuit  where  an  electrical  resistor  must  be  used  to  attempt  to 
reproduce  the  water  effect.  A  combined  resistance  and  capacitance  balancing 
circuit  is  provided  within  the  bridge  circuit  to  assist  that  available  in 
the  amplifier. 

As  mentioned  previously,  good  linear  response  has  been  established 
for  this  probe  in  operation  over  the  range  of  +  0.20  radians  ample  for 
the  usual  model  wave  system  and  relatively  good  insensitivity  with  changes 
in  wave  height  over  the  range  of  *  l*j  inches  required  for  the  present  tests. 

Care  is  taken  to  keep  the  wires  clean  and  at  constant  tension.  No  difficulty 
with  this  arrangement  has  been  njted  due  to  variations  in  surface  film,  nor 
does  the  exact  spacing  of  the  wires  seem  to  be  a  critical  factor  as  has 
been  noted  by  others11.  This  might  be  a  fortuitous  result  of  the  relatively 
high  exciting  frequency  of  2KC  provided  by  the  Brush  amplifier. 

The  foregoing  statements  refer  to  the  characteristics  noted  under 

static  conditions  and  do  not  include  any  effects  due  to  dynamic  action,  such 

as  meniscus  oscillation  at  the  frequency  of  the  waves  being  measured.  Although 
13 

Pearlman  ran  a  series  of  such  tests  on  a  similar  device  and  found  no  important 
effects  at  the  frequencies  of  importance,  the  present  device  should  be  sub¬ 
jected  to  a  similar  series  of  tests  for  the  sake  of  completeness. 
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APPENDIX  B:  COMPUTER  PROGRAM 


The  computer  program  for  evaluation  of  the  wave  resistance  and 

corresponding  wave  spectrum  from  the  lateral  wave-slope  record  was 

originally  written  for  an  IBM-1130  computer  and  has  since  been  modified 

to  run  on  a  GE  235  time-sharing  computer.  This  program  is  shown  in 

Fig.  B-l,  and  is  in  the  FORTRAN  language.  The  basic  equations  to  be 

computerized  are  (18)  and  (19). 

The  data  is  composed  of  approximately  250  wave-slope  values 

(ZETAY)  which  are  read  from  the  recorder  tape  (see  Fig.  3)  into  storage, 

along  with  values  for  waterline  length  (LWL),  ship  speed  (VEL),  the  time 

increment  between  data  points  (DELT),  wetted  surface  (WETSUR),  and  the 

number  of  data  points  (NOZETA).  The  program  first  calculates  and  stores 

values  of  W(NU)  for  fifty-one  values  of  U(NU)  spaced  equally  from  0  to  5 

in  steps  of  0.10.  The  latter  is  controlled  by  choice  of  the  parameter 

k  b  (CAYOB)  which  is  related  by  Equation  (23)  to  an  effective  tank  width 
o 

b  (WID).  The  following  computations  are  then  made: 
ef  f 

NO  2.  ETA 


SY(nu)  =  ^ZETAY  (J)*S1N  (W(^U)*J*DELX) 


l 


(B-l) 


NO  2-ETA. 

CY1nu)=  yZETAVil)*COS 


where 


T=  l 

PELX  =  DELT*VEL 


Successive  values  of  the  sines  and  cosines  are  computed  using 
trigonometric  identities,  instead  of  always  depending  on  the  subroutines, 
to  save  computational  time. 
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1000  SFILE  TESTl,  TEST2.  TEST3 
1100  DIMENSION  U<  100)  ,  V<  100) 

1200  READ,  C/70B,  CALI  B 
1300  DIMENSION  ZETAY<250) 

1409  READ,  DELT,NOZ ETA 
1500  PRI NT”WHI  CH  SHIP  TEST'* 

1600  10*1 NPUT,M 

1700  READ(M)  VEL,  LVL, WETSUR 

1800  READCM)  ZETAY 

1900  CALI  B=8  5»  0 

2000  CAY0=32. 16/VEL**2 

2100  WID=CAY0B/CAY0 

2200  FN=VEL/SQRT(32.2*LWL) 

2300  VK*VEL*3600./6080. 

2400  PRI NT”SHI P  SPEED  IS  ",VK,"  KNOTS," 

2500  PRI NT"FROUDE  NO.  I S", FN 

2600  PRI  NT"WI  DTH  I  S",  WI  D 

2700  N0NU=0 

2800  DO  100  N= 1 , 5 1 

2900  NONU=NONU+ 1 

3000  XNU=N- 1 

3100  U(N)=3» 141 59*XNU/CAY0B 
3200  W(N)  =  (  . 5+C .25+UCN)**2)** . 5) ** . 5 
3300  100* CONTINUE 
3400  RW=0 • 

3500  DELX=DELT* VEL 

3600  C0N=20C9«*DELX**2/( WETSUR* CAY OB* CAL I B**2) 

3700  PRI NT”CON  IS”, CON 

3800  PRINT,  »t"NU  W(NU)  U(NU)  AMP  THETA  PHASE  WEIGHT 
3900  DO  400  NU= 1 , NONU 

4000  sy=o.;cy=o. 

4100  SINDEL=SINCDELX*W(NU)*CAYO) 

4200  COSDEL=COS(  DELX*  W<  NU)  *CAYO) 

4300  S 1  NX = 0 • 

4400  CGSX=1. 

4500  DO  300  J= 1 » NOZETA 

4600  CY=CY+ ZETAY C J) *COSX 

4700  SY=SY+ ZETAY (  J)  *SI  NX 

4800  COSI  NE=  COSX*COSDEL~  SI  NX* SI  NDEL 

4900  SI NE=SI  NX*COSDEL+COSX*SI NDEL 

5000  COSX=COSI NE 

5100  300  *  SI NX=SI NE 

5200  THETA=  180 . *ATAN( U< NU)  /W(NU)  )  /3. 14159 
5300  PHASE= 180 • * ATAN( SY/CY)/3. 141 59 
5400  AY  =  SY  *  SY  +  CY*  CY 
550C  AMP=SQRT( AY) 

5600  NUU=NU- 1 

5700  WEIGHT=1./(<W(NU)**4)*(2.*W<NU)**2- 1.)) 

5800  R=CON*AY*WEI  GHT 
5900  RW=RW+R 

6000  PRINT  1,NUU,  W<NU),U<  NU),  AMP,  THETA,  PHASE,  WEI  GHT,  RW 
6100  1*F0RMAT(I2,2F6.3,3F6.2,F6.4,F1L. 7) 

6L'00  400:  CONTINUE 

6300  PRI NT"SHI P  RESI STANCE  IS",RW 

6400  GO  TO  10 

6500  STOP 


Figure  B-1 :  Computer  Program 


RESI STANCE" 
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CY(No)and  SY(N\>)are  the  components  of  the  spectrum  which  is 

printed  out  in  terms  of  the  magnitude  (AMP)  and  the  phase.  These 

results  are  used  to  calculate  the  wave  resistance: 

NCJNU 


2  (DELX) 


[syM+ 

■[aw 

0 

(W(nu 

>M 

\ ^ 

1 

L-0 

(B-2  ) 


NU  =  0 

These  are  then  printed  out  as  individual  and  cumulative  values  for  each  NU. 

The  total  running  time  for  this  program  on  the  GE-235  is  23  seconds 
per  resistance  point  plus  13  seconds  for  compiling  or  a  total  of  36  seconds, 
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